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Abstract: We consider an extension of the Standard Model with an arbitrary number 
of heavy quarks having general couplings to the Higgs boson. We construct an effective 
Lagrangian integrating out quarks that are heavier than half the mass of the Higgs boson 
and compute the Wilson coefficient for the effective gluon-Higgs vertex through NNLO. We 
apply our result to a composite Higgs model with vector-like quarks coupling to the third 
generation quarks. In the heavy quark-mass approximation, we show that the suppression 
of the leading-order cross section with respect to the Standard Model does not depend 
on the number of vector-like multiplets introduced. We analyse the effects of QCD and 
electroweak corrections through three loops, as well as bottom-quark contributions through 
two loops. 



1. Introduction 



The discovery of the sector responsible for electroweak symmetry breaking (EWSB) is 
the main goal of the LHC. In the minimal description provided by the Standard Model 
(SM), the electroweak symmetry is broken by an SU(2)l doublet that acquires a vacuum 
expectation value. After EWSB, three of the four degrees of freedom associtated to the 
Higgs doublet provide the longitudinal modes of the electroweak gauge bosons, and only one 
real degree of freedom survives: the Higgs boson. In this minimal description, the couplings 
of the Higgs boson to the SM particles are fixed by the mass of the particles themselves. 
Indirect experimental bounds point towards a relatively light Higgs, with a mass of a few 
hundred GeV. On the other hand, a fundamental scalar receives quadratically divergent 
contributions from radiative corrections. A light, fundamental Higgs boson then requires 
very fine-tuned cancellations between the tree-level mass and the higher-order corrections 
to it. 

New physics scenarios try to address this problem introducing some mechanism to 
protect the mass of the Higgs boson. One of the most studied examples is super symmetry. 
In supersymmetry new particles are introduced, in such a way that their contribution to 
the self-energy of the Higgs boson counterbalances the contribution from the SM particles. 
Another possibility is that the Higgs boson is not a fundamental scalar, but a composite 
state of some new, strongly interacting sector [1,2]. Compositeness would explain the 
insensitivity of the Higgs boson mass on ultraviolet physics, as loop contributions are 
cut-off at the compositeness scale. Such a scale is expected to be of the order of a few 
TeV. Higher values would reintroduce fine-tuning problems, while a too low scale would 
be hard to reconcile with electroweak and flavour physics constraints. In this framework, 
the mass of the Higgs boson can be kept naturally light by identifying the Higgs with the 
pseudo-Goldstone boson of some spontaneously broken global symmetry of the new sector, 
in analogy to what happens for the pions in QCD. Quarks masses can be generated by 
mixing the massless fundamental quarks of the SM with heavy composite states of the 
new sector. This mixing induces a coupling of the SM fermions to the Higgs boson. As 
in the SM, only heavy quarks couple significantly to the Higgs boson, as they are mainly 
composite. However, the strength of this coupling is not simply proportional to the mass, 
but it depends on the details of the model. 

The introduction of new particles and the modification of the Higgs couplings can 
change in a significant way the Higgs-boson phenomenology. At hadron colliders, the 
main mechanism for the production of the Higgs boson is gluon fusion. This process is 
mediated by heavy-particle loops and therefore affected both by the introduction of the 
new quarks, and by the modification of the Higgs couplings. Within the SM, the gluon- 
fusion channel has been studied thoroughly. The inclusion of next-to-leading order (NLO) 
[3, 4] and next-to-next-to-leading order (NNLO) [5-7] corrections was necessary in order 
to match the accuracy expected by the experiments and to achieve a good converge in the 
perturbative expansion of the cross section in the strong coupling constant a s . At the LHC, 
these contributions yield an increase of a factor of 2 in the Higgs boson production cross 
section. In composite Higgs models, the LO production cross section is expected to be fairly 
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suppressed with respect to the SM value [8,9]. The actual suppression factor depends on 
the details of the model. For an SO (5) — > SO {A) symmetry breaking pattern in the strong 
sector, with composite fermions embedded in the fundamental representation of 50(5) and 
a global symmetry breaking scale / = 500 GeV, the Higgs boson production cross section 
is expected to be 35% of the SM value. A more detailed study of the Higgs production 
cross section and branching ratios in the different channels was carried out in [10]. It is 
interesting to analyse the effect of higher order corrections in this class of models, both to 
check if their effect is the same as in the SM and to reach the same accuracy as the SM 
predictions. 

In this paper, we construct an effective Lagrangian integrating out the heavy quarks, 
for which we assume a generic coupling to the Higgs boson. We compute the Wilson 
coefficient for the gluon-Higgs effective interaction through NNLO. We apply our result 
to the composite Higgs model described above, where we introduce one or two multiplets 
of heavy top-partners. We include in our study the full bottom-mass dependence through 
NLO [4,11], the two-loop electroweak corrections [12] and the corresponding three-loop 
mixed QCD and electroweak corrections [13]. These effects are implemented in the program 
iHixs [14]. We combine them with the NNLO Wilson coefficient that we compute in order 
to obtain the Higgs production cross section through NNLO in composite Higgs models. 



2. The effective Lagrangian 

We extend the SM quark sector through new heavy quarks that transform under the fun- 
damental representation of the QCD colour group. The number of heavy quarks, including 
the top, is rih- The number of light quarks is n\. We take the light quarks to be massless 
and not to couple to the Higgs boson. We assume instead an arbitrary coupling Y\ of the 
heavy quarks to the Higgs boson. In the SM, Y{ = rrii/v, where v « 246 GeV is the vacuum 
expectation value of the Higgs boson and rrii is the mass of the quark. The Lagrangian 
describing this model is 

"ft n h 

C = C^ CD + Y,^i(^-mi)^ + C Y , with C Y = -Hj2 Y i^i- (2- 1 ) 

i=l i=l 

Here is the covariant derivative in the fundamental representation of the colour group 
and £q CD is the QCD Lagrangian with only the n; flavours of light quarks. 

We focus on the Higgs production from gluon fusion mediated by loops of heavy quarks. 
When the particles that couple to the Higgs boson are heavier than half the Higgs boson 
mass, we can integrate them out. In this limit, we can replace the original Lagrangian (2.1) 
with an effective Lagrangian 

C £ ff = C e ^ 1 — C\H Oi . (2.2) 
The dimension-four local operator 0\ reads [15] 

Oi = \ G'%G'^ , (2.3) 
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where G"^ v is the field strength tensor in the effective theory, and C\ is the corresponding 
Wilson coefficient [16]. The effective Lagrangian £>qq£, 1 appearing in Eq. (2.2) describes the 
interactions among the light degrees of freedom. It has the same form as CJq CD , but with 
different coupling and field normalizations to account for the missing contributions from 
heavy quarks loops. The parameters in the effective theory are related to the parameters 
in the full theory through multiplicative decoupling constants Q. The derivation of the 
decoupling constants is reviewed in [17,18]. 



3. Details of the calculation 



We compute the Wilson coefficient C\ up to three loops. We start from the bare amplitude 
M^gg^H f° r t ne process gg — > H in the full theory, 

M'gg^H = ^(PLK^C • (3-1) 

Here p\ and P2 are the momenta of the two gluons and e^, are their polarizations. 
This amplitude is related to the bare Wilson coefficient C® by [17] 

C 3 Cl - {N 2_ 1){d _ 2){pi . p2) 2 (P1,P2)| W=W M) » ( 3 - 2 ) 

with N the number of colours and d = 4 — 2e the dimension of space-time. The superscript 
"0" denotes bare quantities. The factor £3 is the bare decoupling coefficient that relates 
the bare gluon field G'°' a in the effective theory to the bare gluon field G^ a in the full 
theory, 

y0,a _ /x) ^,0,a 



GT = y/ (3-3) 

We generate the Feynman diagrams T for the amplitude through three loops using 
QGRAF [19]. Diagrams containing two different heavy mass scales appear for the first 
time at the three-loop order. We then expand all diagrams in the external momenta j»i,P2 
by applying the following differential operator [20] to their integrand: 

00 

VT = J2(Pi-P2) n PnT] Pl=P2=0 , (3.4) 

n=0 

with 

A. = i, f, = in,, v, = - m _J d(d - 2) {□,.□, -«tf u }, (3.5) 

and Djj = dp T dpV ■ Differential operators of higher orders are not needed for the expan- 
sion in the external momenta at leading order. 

After Taylor expansion, all the Feynman diagrams are expressed in terms of one-, two- 
and three-loop vacuum bubbles. We reduce these integrals to a set of five master integrals 
using the algorithm of Laporta [21] and the program AIR [22]. The topologies and the 
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master integrals that we find are the same as in the calculation of the Wilson coefficient 
for an arbitrary number of heavy quarks with top-like Yukawa interactions of Ref. [18], 



h = 
h = 



d d h 1 
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(3.11) 



Vi 

v 2 
v 3 

V 4 = (h -k 2 + k 3 ) 2 

r 5 = (h - k 2 ) 2 , 
v 6 = (k 2 - k 3 ) 2 , 

The last two master integrals contain two heavy quarks of different mass, rrii and rrij. 
They have been computed in [23] . The one-scale master integrals can be found for example 
in [24]. 

For completeness, we report the Wilson coefficient in terms of the bare parameters in 
the full theory, 
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Here we introduced the notation 



with 



S t = e~^ E (4vr) e , 

rih n h 

L? = J>g(m°) , L° = J>g 2 K) , 

i=l i=l 

T8 = 5>?, T? = 5>?log(m?), T° = ^ y °log 2 K), (3.13) 

i=l i=l i=l 



Y° 

= • (3-14) 



This expression reproduces the result of Ref. [18] when Y® = 
3.1 Decoupling and renormalization 

The RHS of Eq. (3.2) is expressed in terms of the bare coupling constant a° in the full 
theory and of the bare masses and Yukawa couplings of the heavy quarks, mj and Yj°; 
Ci = Ci(a®, m-, ). The bare strong coupling in the full theory is related to the bare 
strong coupling in the effective theory a' s ° by the decoupling constant (® [17,25], 

«;° = (C 9 °) 2 a° • (3.15) 

The decoupling constants of the strong coupling and of the gluon field (Eq. (3.3)) in the 
presence of an arbitrary number of heavy quarks have been derived in [18], 
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After decoupling, the bare Wilson coefficient is a function of the bare parameters in 
the effective theory and of the bare masses and Yukawa couplings of the heavy quarks in 
the full theory, C° = C°(a'°,m° , Yp). The bare parameters are related to the renormalized 
ones through multiplicative renormalization constants Z{ as 



a' s ° = ^Z' a a'M 



my 



(3.17) 
(3.18) 



The primed parameters in Eq. (3.17) are in the effective theory and the unprimed pa- 
rameters in Eqs. (3.18) are in the full theory. In the MS scheme the strong coupling 
renormalization constant is 



Z' 



1 



+ 



7T 



7T 



2e 



(3.19) 



where f3' , (3[ denote the first two coefficients of the j3 function in the light-flavour theory, 
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The mass renormalization constant reads [26] 
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(3.21) 



where rif is the number of active flavours. In the effective theory rif = ri\ and in the full 
theory rif = ni + rih- Finally, the Yukawa coupling Y® renormalizes as the mass, 



Z Yl = Z mi . 

We review here the proof of this relation, following Ref. [27]. 
Let us consider the bare scalar current 



(3.22) 



(3.23) 



We relate it to the renormalized scalar current j(x) through a multiplicative renormaliza- 
tion constant Zj, 

f(x) = Zjj(x) . (3.24) 
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Let us recall that the bare and renormalizeel two-point Green functions for the fermions 
with an insertion of the scalar current are defined respectively as 



G ' (2) {pi,p 2 ) = i 2 Jdx 1 dx 2 e i( P lxl -P^(T^{x 1 )j (0)iP-(x 2 )) 



(3.25) 

G {2) (Pi,P2) = i 2 Jdx 1 dx 2 e l ^~ P2X ^(T^ l (x 1 )j(0)^ t (x 2 )) 

= i 2 ^dx x dx 2 e i{ ^- p ^ ^J[V^}e i f dxC j{0)^ i (x 1 )^ i (x 2 )^ . (3.26) 

Here C°, C are respectively the bare and renormalized QCD Lagrangian and f[V§°], f[V&] 
denote the functional integral over the bare and renormalized fields appearing in them. 
From the second line of both these equalities we see that G°'( 2 \ are related to the 
bare and renormalized quark propagators 



as 



S°(p) = i jd A xe tpx (T^(x)^(0)) = i jdxe ipx [V^ }e^ dyC °^(x^-(0) 
S{p) = i Jd 4 xe ipx (Ti; t (x)lp t (0)) = * jdxe ipx [V<5>]e l $ dy c ^(x)^^ , (3.27) 

G o,(2) {pjp) = _J_ s o {ph G m^ p) = _A_ S {p). (3.28) 



Therefore 



£0,(2) = p_ G (2) (3>29) 



where Z2 is the renormalization constant of the quark wavefunction, 

V>° = • (3-30) 

On the other hand, comparing Eqs. (3.25) and (3.26) we have 

G°- (2) = Z 2 ZjG {2) , (3.31) 

and therefore 

Zj = Z-) . (3.32) 

The bare Yukawa Lagrangian C Y has the same form of Ly (Eq. (2.1)), but it contains bare 
fields and Yukawa couplings. Treating the Higgs field as an external static field, we can 
use the result for the renormalization of the scalar current j°(x) to relate the bare and 
renormalized Yukawa Lagrangians. We obtain 

Z Yi = Zj 1 = Z mi ■ (3-33) 



-7- 



We finally renormalize the bare Wilson coefficient Ci(a' s ,m,i,Yi) using [15,28,29] 
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Our final result for the renormalized Wilson coefficient reads 
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Similarly to Eq. (3.13), we defined 
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(3.37) 



For Yi = — we recover the Wilson coefficient for an arbitrary number of heavy quarks 
with Standard Model-like Yukawa interactions of [18]. 
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4. Composite Higgs models 



In this Section we review briefly the composite Higgs model of Refs. [30,31]. 
4.1 The Higgs sector 

At low energy, the strong sector responsible for EWSB in composite Higgs scenarios can 
be described by a non-linear sigma model. Such description allows to decouple the light 
pseudo-Goldstone boson from other heavy composites of the new sector. We will consider 
the global symmetry breaking pattern SO (5) — > SO (4), as this is the minimal pattern 
that includes custodial symmetry. An additional U(l)x symmetry is required in order 
to generate the correct Weinberg angle. The SM electroweak group SU(2)l x U(1)y is 
embedded into 50(4) x U(l) x ~ SU(2) L x SU(2) R x U(l) x so that the hypercharge is 
given by Y = + X [33, 34]. We denote by / the scale at which the global symmetry 
is broken. This scale is assumed to be larger than the EWSB scale v ~ 246 GeV. Too 
large values of / would introduce a substantial fine-tuning of the model [32] ; on the other 
hand, if the scale of new physics is too low, large contributions to electroweak parameters 
and flavour physics are expected. We set / = 500 GeV, which corresponds to a ~ 10% 
fine-tuning [32]. 

The effective theory becomes strongly coupled at a scale ~ 47r/. Above this energy, 
a more fundamental description needs to be introduced. If the coupling constant of the 
strong sector g p is not maximal, i.e. g p < 4ir, hadronic bound states appear below the 
strong-coupling scale. The typical mass of these resonances is m p ~ g p f and acts as a 
lower, effective cut-off for the effective-theory description. Thorough this work we will 
assume the presence of such weakly-coupled new states and take as cut-off A = m p ~ 2nf. 

The 50(5) — > 50(4) breaking is realized through a scalar field E(x) subject to the 
constraint 

S 2 (x) = 1 . (4.1) 

In the non-linear representation 

E(x) = £ e n(a,)// , n(x) = -iT & h & (x)V2 , (4.2) 

where So = (0,0,0,0, 1) is the vacuum state that preserves 50(4), T a are the four broken 
generators and h a the corresponding Goldstone bosons. Eq. (4.2) can be rewritten as 

£ = ( ^ sin ^,cos ^] = (E, S 5 ) , (4.3) 



h r f 

where h = ((j) * , — (f)~ , (f) + , 4>°) transforms under the fundamental representation of 50(4) 
and h = \J (/i a ) 2 . The SM Higgs doublets with hypercharge +1/2 and -1/2 are respectively 
<& = (</> + ,(/>°) T and 4> = ia 2 &*. The requirement 

m w = ~j- sm 2 J = (4.4) 
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for the mass of the W boson yields 

Higgs compositeness and the requirement for canonical normalization of the kinetic 
term lead to a rescaling of the physical Higgs field by a factor c a = y/l — v 2 / f 2 . This 
rescaling yields a reduction by a factor c a of the couplings of the Higgs with the electroweak 
gauge bosons. The coupling of the Higgs boson to the fermions is also suppressed, but the 
suppression depends on how the fermions are embedded into representations of 50(5). We 
will discuss this aspect in Section 4.2. 

The reduction of the couplings between the Higgs and the gauge bosons leads to some 
sensitivity of the Peskin-Takeuchi 5 and T parameters [35,36] on the cut-off A. In the SM 
the Higgs boson regulates the logarithmic divergencies of the gauge bosons self-energies. 
The reduction of the Higgs boson couplings to the gauge boson spoils this effect. The 
result is a positive contribution to the 5 parameter and a negative contribution to the T 
parameter [32]. The strong dynamics can further affect electroweak precision observables 
through some higher-order operator. Custodial symmetry is included to protect T, while 
the 5 parameter receives a further positive contribution. We refer to [32] for a more 
complete discussion of these effects. 

The shifts to the 5 and T parameters from Higgs compositeness and from UV physics 
make the model incompatible with EWPT [30-32]. However, other composite states might 
lie below the cut-off of the effective theory and contribute to these observables. In par- 
ticular, quark masses arise through mixing of the SM elementary fermions with fermionic 
bound states of the strong sector. We analyse this scenario and the effects on the elec- 
troweak parameters in the next Section. 

4.2 The fermionic sector 

We extend the 50(5) symmetry of the strong sector to the fermion sector and assume 
that composite top-partners lie below the cut-off of the effective theory description. Many 
viable choices for the embedding of these states into representations of 50(5) have been 
considered, both in the context of five-dimensional and effective four-dimensional mod- 
els [8,32-34]. Following [31], we include fermionic multiplets ^ that transform under the 
fundamental representation of 50(5). Their SO (5)-invariant mass Lagrangian is 

-£ so{5) = rr^¥¥ + /^/OrSXEtv^) , (4.6) 

where the indices i,j allow for the possibility of more than one set of fermionic composites 
and the brackets in the second term indicate the contraction of the 50(5) indices. In 
Eq. (4.6), hermitian matrix and £ is the Higgs field defined in (4.3). 

For models with one multiplet, only a small region of parameter space is allowed by 
electroweak precision measurements [30,31]. The introduction of more sets of fermionic 
composites is inspired by the five-dimensional model presented in [37], which does not 
seem to suffer from such severe constraints. In [31], it was shown that a four-dimensional 
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composite Higgs model with two multiplets of fermionic resonances is compatible with 
EWPT in large regions of parameter space. 

A vector (5) of SO(5) decomposes under SO(4) ~ SU(2) L x SU(2) R into a bidoublet 
(Q,X) plus a singlet T, 



* = (Q,X,T) (5) = (4) (1) ~ (2,2) (1,1) . 



(4.7) 



We assign to *!> an U(l)x charge of 2/3. In this way the SU(2)l doublets Q and X 
get hypercharge 1/6 and 7/6, respectively, and the singlet T acquires hypercharge 2/3. 
Therefore the SM quarks qi and t R have the same quantum numbers as Q and T. We can 
write for them an interaction Lagrangian 



-dm = m\q L Q R + m R T\t R + h.c. 



(4.8) 



The doublet X introduces a new quark of electromagnetic charge 2/3, which mixes with 
the top after electroweak symmetry breaking, and an exotic quark of charge 5/3 that does 
not couple to the Higgs boson. 

Using Eqs. (4.6) and (4.8) and expanding cfp in (4.3) around its vacuum expectation 
value, = ^j^, we obtain the mass terms and the Yukawa couplings of the quarks. The 
mass matrix for the quarks of charge 2/3 reads 
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(4.9) 



The indices u and d denote the upper and lower components of a doublet, respectively. In 
the case of more fermionic resonances, the mass matrix is to be understood in block form. 

Since the mass of the bottom quark is small, we do not expect large effects from 
bottom compositeness. Instead of generating a bottom mass introducing additional £0(5) 
multiplets, we adopt a minimal description and introduce an explicit 5'0(5)-breaking term 



C b = \ b q L ^b R 



(4.10) 



to give a mass to the bottom quark. Therefore the mass matrix for the quarks of charge 
-1/3 is 
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The Yukawa couplings of the top-like quarks are 
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(4.12) 
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where <8> denotes the matrix tensor product. 

The most important bounds on the model come from the S and T parameters and 
the anomalous ZbjJiL coupling. We use the fit to these three quantities that was employed 
in [31]. The scan over the parameter space is done as in Ref. [47]. We choose / = 500 GeV, 
m,H = 120 GeV, and we set the top and bottom masses to [38, 39] 

m t = 172 GeV and m b = 4.16 GeV . (4.13) 

Direct experimental searches impose lower limits on the masses of the new quarks. 
These analyses however assume that the new quarks decay entirely through one specific 
channel {V -> tW~ [40,41], b' -> bZ [42], i! -»• bW+ [43], f -»• qW+ [44], X -»• tW+ [45]). 
Studies carried out in the context of a four-generation Standard Model show that the 
bounds can be significantly lowered when multiple decay channels are open [46]. Follow- 
ing [47], we impose the limits 

m 5/3 > 365 GeV , m 2/3 , m_ 1/3 > 260 GeV (4.14) 

on the masses of the new quarks of charge 5/3, 2/3 and —1/3. 



5. Higgs production in composite Higgs models 
5.1 General LO results 

We first compute the contribution from the charge 2/3 quarks to the LO Higgs production 
in the heavy-mass approximation (m q > 2m#). This is an interesting analysis, as we can 
prove that in this approximation the cross section is suppressed with respect to the SM 
result by a factor that only depends on the scales of the electroweak symmetry breaking v 
and of the global symmetry breaking /. Such result is already known for the case of one 
multiplet [8]. We show that it holds for any number of multiplets. 

Denoting by u^^ SM ^ the LO production cross section in the heavy quark approxima- 
tion in the composite Higgs (Standard) Model, the suppression factor reads 



9 n SM 
u app 



cos(2(fe)//) l 2 _ (1-24)' 
cos ((/>)//) J 1-4 



(5.1) 



For our choice of parameters, R g = 35%. 

Our proof follows the one of Ref. [8]. In the heavy-mass approximation, the SM Higgs 
production amplitude M.g^ H can be written as (Eqs. (3.1), (3.35)) 

^H = f(e, P )^ = f(e, P ) 1 -. (5.2) 

•"'top <J 

In this expression, f(e,p) contains the dependence on the polarization and momentum of 
the external gluons and m top ,Y top are the mass and Yukawa coupling of the top quark, 
respectively. In the presence of more heavy quarks of mass rrii and Yukawa coupling Yj,, 
this result generalizes to 
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so that 



Here M and Y denote respectively the matrices of masses and Yukawa couplings. Using 

Y-™- (55) 

Y - d(h) ' (5 - 5j 

we can rewrite (5.4) as 

/ 3M \ d d 

fl » /2 = " Tr ( M ~W =v m Tl]ogM = v m' ogaetM - (5 ' 6) 

The dependence of the determinant of the mass matrix on (h) (i.e. on s a , c a ) is of the form 

det M = s a c a £(m L ,m R , m-q,, fi, /) , (5.7) 

where ^(mi, mj?, m^, fi, /) is a function of the parameters indicated. We derive the expres- 
sion for det M and give the explicit form of the function £ in Appendix A. Inserting (5.7) 
into (5.6) and using the definition of s a (Eq. (4.5)), we obtain 



sin 



««(2W//) 

cos ((/»)//) • 1 • j 



As anticipated, this result does not depend on any of the parameter-space details, including 
the number of fermionic multipets. 

Finite-mass corrections, bottom-quark and electroweak effects modify this result al- 
ready at LO. In particular, bottom-quark and electroweak contributions are more signifi- 
cant than in the SM. In the SM, the inclusion of the bottom quark lowers the LO Higgs 
production cross section by about 7%, while electroweak effects give a 5% increase. As we 
have seen, in the composite Higgs model the contribution from heavy quarks is strongly 
suppressed. On the other hand, in our description we couple the bottom quark directly 
to the Higgs boson. As a consequence, its Yukawa coupling is reduced only by about 13% 
with respect to the SM value 1 . One therefore expects a larger reduction of the cross section 
from bottom quark loops, of the order of 10%. Similarly, the couplings of the gauge bosons 
to the Higgs are reduced by a factor c a ~ 87% with respect to their SM value. At LO in 
a s , the contribution from electroweak corrections should therefore yield an increase of the 
cross section by about 7%, against the +5% of the SM. These estimates are confirmed by 
the exact numerical values for the LO production cross section that we report in Table 1. 



1 This suppression is related to the mechanism that we adopt to give a mass to the bottom quark. In 
particular, it can be modified in a scenario where an additional SO(5) multiplet is introduced and the 
bottom quark acquires a mass in a similar way as the top quark. 
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a LO \pb] 






SM 


8.8 


8.1 


-7% 


9.2 


+5% 


8.6 


CH 


(2.9 -7- 3.2) 


(2.6 -T- 2.8) 


-10% 


(3.2 -T- 3.4) 


+7% 


(2.9 -r- 3.1) 



Table 1: Leading-order gluon- fusion cross section in the SM and in the composite Higgs model of 
Section 4 for the 7 TeV LHC. The notation (x m i„ -j- x max ) indicates the range of values that the 
quantity x can assume. We report the cross section due to charge 2/3 quarks only (including 
finite- mass effects) , and analyse how it changes with the inclusion of bottom-quark and electroweak 
corrections (<r^ b ° and <t£,°). In the last column we give the total LO cross section a LO . 



5.2 Precise prediction through NNLO 

We now compute the Higgs production cross section in the composite Higgs model through 
NNLO. We include the contribution from the heavy quarks retaining the full mass de- 
pendence through two loops. The NNLO corrections are computed in the effective theory 
approximation according to the Wilson coefficient (3.35). Since all the heavy quarks are 
integrated out from the low-energy effective theory, the only difference with the SM calcu- 
lation is in the expression of the Wilson coefficient. The remaining part of the calculation 
is the same as in the SM [5-7]. Following the approach of [13], the NNLO corrections are 
normalized to the exact LO cross section according to 

(„N 'N LO;heavy \ 
° 1 fee cA 

-LO-Mavy J " 
/ effective 

Since bottom-quark effects are more important than in the SM, we compute them exactly 
through NLO [4, 11]. We also include the full two-loop SM electroweak corrections of 
Ref. [12] and the three-loop mixed QCD and electroweak corrections derived within an 
effective-theory approach in Ref. [13]. Both the electroweak corrections are rescaled by 
the factor c a that reduces the coupling of the Higgs to the gauge bosons. All the effects 
described here are included in the code iHixs [14], which we use for the calculation of our 
results. 



In Table 2 we present the full NNLO cross section in the SM and in the composite Higgs 
model. The results are similar for the case of one and of two multiplets of composite 
fermions. We use the MSTW2008 NNLO parton distribution functions [48] and set the 
renormalization and factorization scales to fx = fif = fj, r = mif/2. The 35% suppression 
factor R g computed in Section 5.1 is confirmed through NNLO. We estimate the uncertainty 



a SM ^ 


4 + i % 


<4ale % 


a CH \pb] 


4 + i % 


41 % 


R'g 


17.6 


+9 % 


-10 % 


(5.9 -7- 6.4) 


+(6 -7- 12)% 


-(7 - 11)% 


(34 -7- 37) % 



Table 2: Gluon-fusion cross section through NNLO in the SM and in the composite Higgs model, 
with the corresponding scale variation errors. The factor R' is defined as in Eq. (5.1), but for the 
full result through three loops. 
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SM 


CH 


a NLO 
a LO 


+ 75% 


+ (77 -=- 78) % 


a NNLO 


+ 106% 


+ (108 -T- 110) % 





Table 3: NLO and NNLO K- factors in the SM and in the composite Higgs model. 

due to higher order corrections by varying the scale ji in the interval (rnff/4, ma)- This 
scale variation uncertainty is similar to the SM one. 

Finally, in Table 3 we compare the K-factors in the SM and in the composite Higgs 
model. As in the SM, the K-factors are large and the NNLO result is about twice as big 
as the LO cross section. 

6. Conclusions 

We presented the construction of an effective theory for extensions of the Standard Model 
with an arbitrary number of heavy quarks coupling to the Higgs. We assumed a general 
form for the Yukawa couplings of these quarks to the Higgs boson. This situation arises 
for example in the context of composite Higgs model, where the mass of the quarks can 
be explained through the mixing of the fundamental SM particles with heavy composite 
fermions. We computed the Wilson coefficient of the effective Higgs-gluon vertex through 
O(al). We used our result to compute the Higgs production cross section through NNLO 
in a composite Higgs model with an 50(5) — > 50(4) global symmetry breaking pattern 
and one or two multiplets of composite fermions transforming under the fundamental rep- 
resentation of 50(5). We showed that, in the heavy quark-mass approximation, the LO 
production cross section is suppressed with respect to the SM value by a factor that de- 
pends neither on the details of the parameter space nor on the number of multiplets. As 
in the SM, the NNLO result is enhanced with respect to the LO cross section by approx- 
imately a factor of 2. The scale variation errors also behave in a similar way as in the 
SM. We included in our result also the full dependence on the bottom quark mass through 
two loops and the two-loop electroweak and three-loop mixed QCD and electroweak cor- 
rections. Both these effects are enhanced with respect to the SM. As in the SM, they give 
contributions of opposite sign, which cancel. 

In this work we applied our result for the Wilson coefficient to a specific beyond-the 
Standard Model scenario, but the calculation can be extended to any model with additional 
new quarks in the fundamental representation of the colour group. 

Acknowledgements 

We thank Babis Anastasiou and Giuliano Panico for many useful discussions and for their 
comments on the script, and Achilleas Lazopoulos for providing the preliminary version 
of iHixs. We greatly appreciated the hospitality of the ETH theory group during parts of 
this work. This research is supported by the DOE under Grant DE-AC02-98CH10886. 



- 15 - 



A. Analytical form of detM 

We derive here the analytical result for the determinant of the mass matrix of charge 2/3 
quarks, detM, for an arbitrary number of multiplets. Let us recall (Eq. (4.9)) that the 
mass matrix for the charge 2/3 quarks reads 



/ 



M 



m, 







m* + £# 4* 

c 2 £l± ™. T . _L.2/M 



'a 2 





s c i± 



\ 



yrriR s a c a &= s » c a^ rn^+cifiij 



(A.l) 



Successively taking linear combinations of lines/columns of M, we can recast it into the 
form 

\ 



M' 





°a 2 










s^+2c^ 



my 



in 



SqCq V2 




(A.2) 



where 



m 



T 







s 2 ^ 

*a 2 U 



B 



and D 




\ / 

Because of the properties of determinant, 





—my 



(A.3) 



det M = det M' = det D det (A - BD" 1 ^ 
1 



m T L W- 1 m,R det(m^) det(//i) det(W) 



where 



IV2 



(A.4) 



(A.5) 



The quantity in square brackets in Eq. (A.4) corresponds to the function £ introduced in 
Eq. (5.7). 
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